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Abstract 

We consider summations over digamma and polygamma functions, often 
with summands of the form (±l) n tjj(n + p/q)/n r and (±l) n ^( m )(n + p/q)/n r , 
where m, p, q, and r are positive integers. We develop novel general integral 
representations and present explicit examples. Special cases of the sums reduce 
to known linear Euler sums. The sums of interest find application in quantum 
field theory, including evaluation of Feynman amplitudes. 
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Introduction 

The evaluation of Euler sums [3, 4, 5, 6, 7, 9, 11, 15, 16] is not only a beautiful 
subject in classical real and complex analysis, but has application in quantum field 
theory. Specifically, in quantum electrodynamics, one evaluates Feynman amplitudes 
from certain integrals [14, 26]. These often lead to Mellin transforms and Euler sums 
[8, 21, 22, 23, 28]. The problem of closed form evaluation of Euler sums has stimulated 
research in symbolic computation and integer relation detection techniques [3, 28]. 

Among many possible variations (e.g., [3, 24]), various one dimensional Euler sums 
may be written as 

n=l u n=l u 

oo \u(r)]p oo f/jMlP 

n=l 11 n=l 11 

where the generalized harmonic numbers are given by 

^ r) -E^ H n = H$;\ ^^EMt-- ( 2 ) 
j=i J j=i J 

A rich combination of elementary and special functions arises in the evaluation of Eu- 
ler sums. Special functions which typically appear are the polylogarithms, gamma, 
digamma, and polygamma functions, zeta functions, Clausen functions, hypergeo- 
metric functions, tangent integral, and logarithmic-trigonometric integrals [16, 3, 9]. 

This paper considers sums over digamma and polygamma functions at rational 
argument, wherein special cases reduce to one dimensional linear Euler sums. The 
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tight connection between harmonic numbers and the digamma function is clear: H n = 
i]){n + 1) — if)(l) = ip(n + 1) + 7, where 7 is the Euler constant. This extends to 
generalized harmonic numbers and the polygamma function, 

' V> (r-1) (n + l) -^ (r-1) (l)] , 

where ^ (r_1) (l) = (— l) r (r — 1)!C (r) and C is the Riemann zeta function [13, 19, 25, 27], 
so that 

**. = ft9,| + 1) - (-ir'cMcfa), r,„ > 1. 

Another set of numbers that often figures prominently in special cases of our summa- 
tions is 

n , 

7 + V(n + 1/2) = -2 In 2 + 2 £ — — . 

Our results extend many of the previous sums calculated by de Doelder [11]. A 
significant physical motivation is provided by the very recent work of Ogreid and 
Osland [21, 22, 23], and we build upon their discussion. These authors have cal- 
culated one and two dimensional Euler series associated with Feynman diagrams, 
invoking properties of the dilogarithm and trilogarithm functions. Our focus is not 
with polylogarithms [20], but to bring to bear especially further properties of the 
polygamma and Clausen functions. 

Previously, we considered several nonlinear Euler sums, including one closely re- 
lated to 5*1,4,2, X^Li Hn/{n + l) 2 [9]. We discussed the application of logarithmic- 
trigonometric integral, Fourier series, and contour integral methods. We conjecture 
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this sum to have the value (859/24)C(6) + 3C 2 (3) [3], where C(6) = tt 6 /945. 

The sum Si t i iq has a closed form, and a history that seems to go back to Euler 
[24]. The companion alternating sums and A' l l q can be written as a polynomial 
in zeta values when the weight 1 + q is odd, while the alternating sum S' l l q has a 
closed form for any weight [24, 15]. In Ref. [24], Sitaramachandra Rao formulated 
twelve classes of infinite series. With his notation, S^i,, = Hi(q) and = —H 7 (q). 
Among many useful results and consideration of a Ramanajun sum, this author gives 
the closed form for the sum S' l l q = H 4 (q). 

We first consider sums over the digamma function, containing summands with 
(±l) n ip(n + p/q)/n 2 . Such sums are amenable to much simplification and several 
explicit alternative forms, owing to the properties of the Clausen function CI2. This 
function is known at rational multiples of 7r and in fact is expressible in terms of 
the trigamma function at rational argument [12]. We next present similar results for 
summands with (±l) n tp(n + p/q)/n, before proceeding to sums over (±l) n ip(kn)/n p . 
For these more general sums we construct new integral representations and illustrate 
their reduction in a few special cases. We then make an extension to sums over the 
polygamma function. The integral representations presented here too subsume several 
previous results. Before concluding, we offer several variations on digamma and 
polygamma series with both integer and rational arguments and make a suggestion for 
the further study of linear digamma and polygamma series with reciprocal binomial 
coefficients. We have relegated to an Appendix a quick estimation of many of the 
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digamma and poly gamma sums of interest. 

Sums over (±l) n ip(n + p/q)/n 2 

We first present a general result for sums over the digamma function at rational 
argument. For integers q > 1 we have 
Proposition 



oo y q-1 ' 

E E 

n=l r =0 L 



7 + ijj \ n + 



q 



q 



q-1 



2nf 



F + ^C(3) + ,gC, n - 



E^E 



n=l 



r=0 







7 + ip \n-\--\ 


= q 







glngC(2), 
(3a) 



7 2 3 



C (3) + /eV^ + " 

j'=i 



+|lngC(2), 

(36) 



sin /c6* 



A; 2 



In p 



-dp. 



(4) 



2 8q 

where the Clausen integral is given by [1, 12, 18] 

QUO) = - / ln[2sin(t/2)]d* = V 

/- 1 , / p sin 9 \ dp . n r 1 

= / tan" 1 — — = - sin 9 / — 

Jo \l — pcosOJ p Jo p 2 — 2pcos0+l 

For proof, we start with the multiplication formula for the digamma function 
[2, 17], 

I & H — = qip{bq) — q In q. (5) 
r =o V 9/ 

We apply this formula with real b put to integer n, multiply by (±l)"/n 2 on both 
sides, and sum from n — 1 to oo, resulting in 



oo i q-1 ' 

E-^E 



n=l n ' r=0 L 



7 + *0 n + - 



oo -y 

<?E — [7 + "0(^)1 -glngC( 2 )> ( 6a ) 



n=l 



5 



and 



\n q-1 



n=l 



n 



r=0 



7 + ^ \ n + 



n=l 



!i "[7 + ^M] + |lngC(2). (66) 



For Eq. (6b), we used the infinite series for the alternating zeta function, Li n (— 1) = 
(2 1_n — l)C(n), where Li g (z) = Y^=\ zU l nq is the polylogarithm. By applying a 
Corollary of Ref. [21], we obtain the Proposition. 

The case of q — 1 is highly degenerate in Eq. (3). Then both the CI2 and logarithm 
terms vanish, resulting in simply £(3) for (6a) and C(3)/8 for (6b). We work through 
a series of additional special cases, due to the applications interest, in order to make 
connections with the literature, and because the Clausen functions may be much 
simplified. 

For q = 2 in Eq. (3), we obtain from the previous q — 1 case 



n=l 



oo 1 

E 

n=\ 
hi) 

n 2 



7 + *0 ( « + 



1 



7 



C(3)-21n2C(2), 



7 + ^(^+2 



-C(3)-27rG + ln2C(2), 



(7a) 
(76) 



where [18] 



G = C1 2 



7T 



£ 



(-i) r 



'o(2m + l) 



16 



[0'(l/4) - 0'(3/ 4 )] ~ 0.915965594177219015, 



(8) 



is Catalan's constant, and 0' is the trigamma function. Since —■0(1/2) = 7 + 2 In 2, 
Eqs. (7) agree with results obtained by de Doelder [11] by way of logarithmic integrals. 
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Equations (7) also result from the direct combination of Eqs. (B.l) and (B.3) of Ref. 
[21]. 

We can just as well apply the method of proof of Eq. (3) with b — > 2n in Eq. (5) 
before summing over n. This gives, for instance, 



OO 



n=l 



rr 



7 + i= (2n + l - 



14C(3) -47rG-21n2C(2). 



(9a) 



This equation also follows from the direct combination of Eqs. (B.3) and (B.7) of Ref. 
[21]. The corresponding alternating sum makes use of a number of the properties of 
the Clausen function, so that we record it as 
Lemma 



E 

n=l 



7 + i> (2n + X - 



14C(3) + C(2) In 2 + 2ttG - 87rCl 2 (7r/4), (96) 



where [18] 



C1 2 (tt/4) = ^-[^'(1/8) - 2(V2 + 1)tt 2 - 8(2^ - 1)G], 



(10) 



In obtaining Eq. (9b), the expression Cl 2 (37r/4) + 2Cl 2 (57r/4) + 3C1 2 (7tt/4) has been 
reduced by using the symmetry Cl 2 (7r + 9) — — Cl 2 (7r — 9) and the duplication 
formula satisfied by the Clausen function. There results Cl 2 (57r/4) = — Cl 2 (37r/4), 
C1 2 (7tt/4) = -C1 2 (tt/4), and C1 2 (3tt/4) = C1 2 (tt/4) - G/2. 
As a follow on to Eqs. (9) we have 



OO 1 

si 



71=1 



rr 



7 + ijj [ 3n + - 



63 

y C(3)-107rCl 2 (7r/3)-21n2C(2), (11) 



where [18] 



C1 2 (tt/3) 



1 



2^/3 



[^(1/3) -2^/3]. 



(12) 



For the corresponding alternating sum we find 
Lemma 



n=l 



7 + tjj 3n + 



1 



= — C(3) +37rCl 2 (7r/3) +ln2C(2) - — ttG - 87rCl 2 (7r/6) 



63 

T C(3) + C(2)ln2-147rG. 



(13) 



In obtaining Eq. (13), relations such as Cl 2 (77r/6) = — Cl 2 (57r/6), Cl 2 (ll7r/6) = 
-Cl 2 (7r/6), (1/2)C1 2 (tt/3) = Cl 2 (7r/6)-Cl 2 (57r/6), C1 2 (tt/6) = (2/3)G+(l/4)Cl 2 (7r/3), 
and Cl 2 (57r/6) = (4/3) G — Cl 2 (7r/6) have been used. In addition, we find 



OO -y 



n=l 



n- 



7 + ip ( An + - 



56C(3) -47rG- 167rCl 2 (7r/4) -21n2C(2), (14) 



and 



~ (-If 

„2 



n=l 



7 + if? (^An + - 



= 56C(3) + 77rCl 2 (7r/4) - 2nG - 167rCl 2 (37r/8) 
-27rCl 2 (77r/8) + C(2) In 2 - 147rCl 2 (7r/8). (15) 



The details of deriving Eqs. (14) and (15) are omitted, although we note that in 
the process we determined that £^°=i { -^f- [7 + ip (An)} = (253/32)C(3) + (n/2)G - 
47rCl 2 (7r/4). 

Similar to the proof of Eqs. (3), we have the following 
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Proposition 



OO 

E ^2 



n=l 



^ + ip [kn + - 



_ 2fe— 1 fc— 1 

-A; 2 C(3)+27r £ iCl 2 (7rjA)-7r^jCl 2 (27riA)-2C(2)ln2, 

(16a) 



E 

n=l 



-l)r 



j'=i 
7 



2fe-l 



7T \ 



-vrEjCl 2 (^ + £)+C(2)ln2, (16&) 

In evaluating such sums, the Clausen function at rational multiples of 7r can always be 
written in terms of the trigamma function at rational argument and the sine function 
at rational multiples of ir. Specifically, we have for integers q > 2 [18], 



ci 2 V 



i q - x 
E 

r=l 



4g2 



smr-7r 



2-7§*'fe) Sin ( r 9"f%^Bi^W2 9 )' 



- 2 ^4 sin(rp7r/g) 



(17) 



(By the use of the discrete sine transform, Eq. (17) can be inverted to yield values of 
the trigamma function in terms of the sine and Clausen functions.) 

Extension to sums over (±l) n ip(n + p/q)/n 



Parallel to the derivation of Eq. (6) from Eq. (5), we have 



9-1 



E-E 



n=l nP r=0 



7 + ip [n + - 

V 9, 



?E —p [7 + ^ Wl - 9 ln ?C(p), Rep>l, (18a) 

n=l n 



and 



-1 



9-1 



E^^E 

n=l r=0 



7 + ijj [n + 



?E 



n=l 



[7 + V(gn)]+glng(l-2 1 " p )C(p), Rep> 1. 

(186) 



When p — > 1, Eq. (18a) diverges, while Eq. (18b) is replaced by 



E^E 



n=l 



r=0 



7 + ip in + 



n=l 



[7 + -0W] + gin gin 2. (19) 



By using a result of Appendix D of Ref. [21], we obtain 



n=l 



r=0 



7 + ijj I n + - 



Cli (^ + ^ 

V ? g, 



where 



Cl 1 (^) = f:^ = -ln|2sin(^/2)|. 



fc=i 



The g = and g = 1 cases of Eq. (20) lead to 



oo (_-t\n 1 

E^[7 + ^H] = ^ln 2 2, 

n=l n Z 



and 



E ( ~ 1] " 



n=l 



n 



- Z C(2) + 2 In 2 2. 



>+g lng 
(20) 

(21) 

(22a) 
(226) 



Generalization to sums over (± 1 ) n ip [kn)/n p 

The generalization to linear Euler sums with higher inverse powers of n is of 
continuing interest to both applications and theory. Accordingly we present new 
integral representations and show that special cases thereof recover known results. 
We have 
Proposition 



E 



r ,/, ,m ±1 r j u% " e " f 1 l ~t K 



dt 



r(p) 



t J (e u T t k ) 



, Re p > 1. 

(23) 
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In the case of the alternating sum, Eq. (23) holds for Re p > 1. Furthermore, when 
k is an integer, this equation may be rewritten as 

^( ±:L )\ ,/, im ±X f°° j uP ~ l ^ f 1 pdt /o^ 

E^ + ^ kn + 1)] = W)Jo (, "T^^L (24) 

With the aid of the functional equation for the digamma function, one has the relation 



2_, nP iT+V^+lj]-^ ^ [7+ ^fc n )j+ < (2 - P _ 1)c(p+1) " bot t om " S i g n 



(25) 

With the change of variable v = t h , the integral over t in Eq. (24) can be written in 
terms of the hypergeometric function [2], 

f 1 Vdt e~ u ( l+ 3 l+j + k \ 
Then changing variable to x — exp(— u), another form of Eq. (24) is 



^(±1)- ±l(-iy- 1 k ^ 1 fUn^x ( l+j l+j + k \ 

(27) 



n=l ^ \r) j=Q 



r(p) ^ J 



Equations (23) and (24) show the special character of the sum when k — 1. In this 
case, the integral over t is elementary, and a change of variable gives the representation 



00 (+l) n (-l) p r 1 ln p ~S; 



(-1) P /"! (I 1 



/ f- + ) ln^ 1 v ln(l - v)dv, "top" sign 

JO \D 1-1)/ 



T(p) Jo \v 1 — v 

m 



■ / ln p_1 u ln(l + u)du, "bottom" sign. 

jo \i> 1 + vJ 
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The subcase at n = 3 recovers results equivalent to those of de Doelder [11] and 
Sitaramachandra Rao [24]: 

oo i f- 

E-i[7 + ^+l)] = -C(4), (29) 

n=l U * 

E L ^-b + 1>(n+l)] = - T C(4) + 2Li 4 (l/2) + i C(3)ln2--ln 2 2 + -ln 4 2, (30) 
with Eq. (30) exhibiting a needed correction to Eq. (8) of Ref. [11]. We also obtain 

OO 1 

E-l[7 + ^(n + l)] = 3C(5)-C(2)C(3). (31) 

n=l n 

Equation (30) shows the needed correction in de Doelder's Eq. (8) in both the 
(7r 2 /12)ln 2 2 and — ll7r 4 /360 terms. The corresponding corrected integral is given 
by 

/; -'"'^-U = 2Li 4 (l/2) - £ + I«3) b2 - ^ 2 - I .„« 2. (32) 

For the case of p an odd integer in the alternating sum of Eq. (28), Flajolet and 
Salvy [15] introduced the constants /i q , with Eq. (30) giving —/ii. Here the explicit 
form of these constants emerges from our novel integral representations (23) and (28). 

For the proof of Eq. (23), we first use an integral representation of the digamma 
function, giving 




We then introduce the representation n~ p = [r(p)] _1 J °° e~ nu u v ~ x du and interchange 
the summation and the integrations. The needed sum is a simply a geometric series, 
and Eq. (23) follows. 
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Extension to polygamma series 

We now discuss the extension to polygamma series. We have 
Proposition 

°o (-\--\\n _|_1 ,00 ,1 f k Wtdt 

V ^L^)(kn + 1) = ±L u p - l e- u du / - ' I -, Re p + j > 1. 

(34) 

For the alternating sum, the range of convergence extends to Re p + j > 0. With a 
change of variable, we have the alternative representation 

> —ibU'fkn + 1) = ' x / ln p L vdv p— , Rep+7>1. 

,S n p ! T{p) Jo Jo {l-t){lTt k v) , F J 

(35) 

By using the functional equation of the polygamma function [1] , we have the relation 



^ nP =-\ nP & ! (2 (P+J) — l)((p + j + 1) "bottom" sign 



(36) 

Equation (34) follows from 

V l^L^0)(A;n + 1) = - Y / j— — In' tdt, (37) 

h nP V K ] ti ^ Jo (1-0 1 ; 

introducing an integral representation for n~ p , and interchanging summation and 
integration. 

The very special case of j — p — k — 1 gives 

£^<» +1 > = {«s)-$.»2 < 38 > 
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The top line of Eq. (38) recovers a known result [23]. For the alternating sum, we 
used [10] 

"i lnyln(l + y) 



Jo 



y 



-dy = 2S lt2 (-l) - Li 3 (-1) + ln2[Li 2 (-l) - C(2)], (39) 



and S 1:2 (-l) = C(3)/8. 
We also obtain 

from [10] 



E("l)V(n) = ~C(2) = 

n=l 4 8 



1 lnt dt 1 /-i / 1 1 \ , , 1 

. n = - / 1 In t dt = - 

lo 1 -t 2 2 Jo VI -f 1 + t/ 2 

More generally, we have for integer j 



f 

Jo 



E(-i) ? 



n=l 



1 \n J t dt 1 r 1 ( 1 1 , , . 

+ ) ln J i dt 



and 



o 1 - t 2 2 Jo VI - 1 1 + t 
_ 1) , [2 -0' + i)_ 1]j!C(j + 1)? 



1 ln i tln(l^t) _ f 1 ln j (l =f w)lnw 



n=l 



y ^±±j-^)( n + 1) = / v ^ ' " y dt = / - v ^ ' ~^"~ dw 

■ 1 Jo I — t 



= 1 

Jo 



w 



= ±l f l ln " 1{lTw) ^ wdw , (43) 
2io 1-ui V ; 

where in the last step we integrated by parts. For Eq. (42) we combined the even 
and odd cases for j of Ref. [17]. We may record a few special-case integrals for use 
in conjunction with Eq. (43): 

/ o 1 ln(l + t)^ = 2[l-C(3)], (44a) 



(40) 



(41) 



(42) 
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f in(i + 1)^^ = r^ 4 - 9 °) ln2 ' ( 446 ) 

Jo 1 — t 15 



and 



r 1 ln 4 f r/f 

/ ln(l + t) — — = -241n2[C(5)-l]. (44c) 
J o 1 — t 

As an extension of Eqs. (34) and (35), we may present a five-parameter summation 



in: 



Proposition 

> —ip (3) (kn+a)z n = / / ln p ^cfc / — 7-, Rep+j > 1. 

^ up r K ' Tip) Jo Jo {l-t){\T zvt k Y F J 

(45) 



On the left side of this equation, the (±l) n factor could just as well be absorbed in 
the new z n factor. The proof of Eq. (45) proceeds as for Eq. (34), again followed by 
a change of variable. The details are omitted. 

Additional digamma and polygamma series 

The many possible extensions of our work include digamma and polygamma series 
with these functions at rational arguments and denominators of the form n a (n + 
l) /3 (n + 2) 7 • • • •. We present in passing a few simple examples: 



and 



t ^f^ 1 /^ = ~IC(S) - 2C(2) + 8 In 2 + 2 C (2) In 2. (48) 

n=l n \ n ' ) Z 
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As an area of possible future research, these sums suggest the investigation of series 
of the form 

±&r b ?«W q)] (49.) 
(") (?) - 

and 

where m, p, q, j, £, ... and a, ft, ... are integers. 
By using the series 

00 t e t 

S ^ + l)« + 2)-« + J) = (J+lji' Fl(1 - 1;J + 2i ' ) - <50) 
and an integral representation for the digamma function, we find 
Proposition for integers j > 1 



f 1 n(n + l)---(n + j) (j + 1)! Jo 



l + --t 2 F 1 {l,l;j + 2;t) 



dt 
l-t' 



(51) 



From the series in Eq. (50), we can multiply by P and integrate from to z to find 

00 z t 1 r z . 

S i(i+l)(£ + 2)...(i + j)(i + j + l) = ^Jo t3fm ' (52) 

where f(t) denotes the right side of Eq. (50). The point here is that we have, 
for instance, [t + (1 — t) ln(l — £)]/£ = (t/2) 2^(1,1; 3; t) so that we can always 
successively obtain the needed hypergeometric functions in Eq. (51) in terms of 
logarithmic functions. 
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Generalizations of Eq. (50) are readily performed. We have 
Proposition for integers j >l 



E 



[7 + ^(fcra + 1)] 



{ n(n + 1) • • • (n + j) (j + 1)! Jo j 



l + --t k 2 F 1 (l,l;j + 2;t fe ) 



dt 
1^? 



and 



E 



^ m \kn + 1) 



1 



V^ l( l,l;i + 2;^ hmftft 



(53) 



(54) 



^in(n + l)-..(n + j) (j + l)!io ' — 1 ' ' i-t 
Just as we derived Eqs. (3) and (18) from the multiplication formula for the digamma 
function, Eqs. (51)-(54) can serve to construct generalizations of Eq. (46) and related 
sums. 

As a set of auxiliary sums complementing the previous sums of this section, we 
present the following 
Proposition for integers j > 1 



E 



= £(2- m -l)C(m+l)+j + l-21n2. 



(55) 



By the use of partial fractions, these sums lead to several related ones. Since these 
are alternating sums with decreasing summand and first term negative, their values 
are negative for all j. 

For the proof of Eq. (55), we first derive 



E 



^ n(n + z) z 
We then differentiate this equation j times to find 



(56) 



E 



(-iy 



\ n(n + z)i +1 



m=0 



ml z i- m+1 



(57) 



17 



We evaluate Eq. (56) at z — 1, separating the m = term of the sum: 



E 



4)n 



^ n(n + 



= E 



m=l 



ml 



2""V ( 



m) 



^ (m) (2) 



+ V 



V>(2). (58) 



To find Eq. (55), we then use the special values ip(2) = 1 — 7, ip(3/2) = 2 — 7 — 2 In 2, 
^Cj)(2) = (_i)i+ij![((j + 1) - 1], and V M ( 3 / 2 ) = (-l) m+1 m!{2 m+1 [C(m + 1) - 1] - 
C(m + 1)}. 

We next indicate an alternative proof of Eq. (55). This method of proof works 
equally well in determining 
Proposition 

-D™ . r j 
— 7 = (-1)' 

m=l 



E 



^nj +1 (n + l) 

In the case of Eq. (55), we put 

00 / 



(-l) m (2' m - l)C(m + 1) + 1 - 2 In 2 



(59) 



I] 



E 



1 



(60) 



^n(n + iy +1 n tl(n+l)i U n + lj' 

Therefore, we easily find the recursion relation Tj = 7}_i + (2~ j — + 1) + 1, with 
T = 1 — 2 In 2. Iteration of this recursion relation gives Eq. (55). 
Similarly, if we put 

-D" ~ (-i) n /1 1 



— = y 

^[rb) +1 (n + l) £r[ rP \n n+1 

we quickly find the recursion relation Sj = —Sj-i + (2~ j — + 1), with S = T 
Iteration of this recursion relation leads to Eq. (59). 



(61) 
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By using the functional equations of the digamma and polygamma functions, Eq. 
(57) can be rewritten as 



E 



E 

m=l 



■1) 



1 



ml zi~ m+l 



2""V (m) 



-m—l 



+ 



:-l) m 2- m ml 



-m—l 



(62) 



Evaluation of this equation at rational values of z offers a multitude of other sums. We 
have ip( n \x) = (— l) n+1 nl((n + 1, x) [1], in terms of the Hurwitz zeta function £(s, a), 
and evaluation of the Hurwitz zeta function at rational argument presents no difficulty. 
For instance, we have V (m) (V 4 ) = (-l) m+1 ml((m + 1, 1/4) for use in evaluating the 
sum J2™=i(— l) n /[n(n + 1/2) J ' +1 ]. By combined use of the reflection and duplication 
formulas for the digamma function, we have simply ^(1/4) = —7 — 3 In 2 — ir/2. 

As a further development, we present 
Proposition For integers N > 1 

1 



OO -y 

§ln?M)('+j) 
(-1)' 1 1 



C(2) + E^( , )[^ + l)+7], 



fc=i 



fc V k 



N 



^ nut +3) ^ c(2) + £ 



l) fc f JV> 
k j 



k " 



^(1+ 2 ) + 



(63a) 



, (636) 



The N = 1 case of Eqs. (63) gives the above Si and Ti sums. The proof of Eqs. (63) 
proceeds from the well known expansion 

-l) fe N 



k I x + k 



We then multiply both sides of this equation by (±l) x /x and sum over x from 1 to 00. 
For Eq. (63b), Eq. (56) is applied, while for Eq. (63a), the widely known summation 



(64) 
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Y^=\[n{n + z)\ 1 = [^(1 + z) + 7]/z is applied. The /c = term is treated by noting 
the expansion [17] 

oo 

^(1 + x) = - 7 + ^(-l)^^)^" 1 - (65) 

fc=2 

This gives, for instance, lim fc ^ ['^'(l + k) + j\/k = C(2), and leads to Eqs. (63). 

The determination of the nonalternating sums corresponding to Eqs. (55) and 
(59) is easily made. We find 



OO 1 

y — - — = f-iy 



m=l 



(66a) 



and 



oo i 3 

These sums may be derived either from a recursion relation or from a summation in 
terms of polygamma functions. 

The sums of this section are intimately linked to hypergeometric summation. For 
instance, we have 3-F 2 (l, 1, 1; 2, 3; t) = 2X)£Li t k ~ l /k 2 (k + 1) and more generally 

oo jk— 1 

m+1 F m (l, 1, . . . , 1; 2, 2, . . . , 2, 3; *) = 2 £ fcm(fc + 1) - (67) 

Of course, the polylogarithm itself is given by Li m (x) = x m+1 F m (l, 1, . . . , 1; 2, . . . , 2; x). 
Particularly in cases of the 3F2 function evaluated at unit argument, the reductions 
of Watson, Whipple, and Dixon [2, 29] can be brought to bear. However, we do not 
pursue this direction here. 



20 



Summary 

Among many developments, we have investigated summations over digamma and 
polygamma functions at rational argument. We have considered both positive and 
alternating series. By exploiting properties of the Clausen function Cl 2 we have pre- 
sented explicit, simplified expressions in many cases. We have developed novel general 
integral representations and given explicit instances. Special cases of the sums reduce 
to known linear Euler sums, or to related digamma and polygamma series. We have 
mentioned connections to the theory of polylogarithms, hypergeometric functions, 
and zeta functions. 

The sums of interest have relevance to calculations in quantum electrodynamics. 
In particular, they support the determination of Feynman amplitudes. Furthermore, 
our sums appear at a sort of cross roads of mathematical physics and analytic num- 
ber theory. Part of the relevance to analytic number theory is the connection of 
subcases to Euler sums, and representations in terms of the Riemann and Hurwitz 
zeta functions. In addition, our results may serve as a type of test bed for symbolic 
computation, and as a stimulus for investigation of further linear sums as in Eqs. 
(49). 
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Appendix: Simple sum estimations 



In the sum 




(Al) 



E- = m(^)+7 + o(l), 



(A2) 



where 7 is the Euler constant. Therefore, for the sum 



n=l 



E 



00 



[y + ip(kn + 1)] 




In £; + 7 + 0(1) 



(A3) 



the leading behaviour is given by the comparison integral f^°(lax)/x p dx = T(2)/(p — 
l) 2 . Therefore there is absolute convergence of the sum of (A.l) for Re p > 1. On 
the other hand, the alternating sum in Eq. (A.l) converges for Re p > 0. 

Similarly, for polygamma series we can make use of the known asymptotic ex- 
pansion [1]. As a very simple example, the leading behaviour of Y,'^=i4> < " m Kkn)/n p is 
given by f£° x~ p ~ m dx. Therefore absolute convergence obtains when Re(p + m) > 1. 
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